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The recently detected polarization of the cosmic microwave background (CMB) holds the potential 
for revealing the physics of inflation and gravitationally mapping the large-scale structure of the 
universe, if so called B-mode signals below 10~ 7 , or tenths of a fjK, can be reliably detected. We 
provide a language for describing systematic effects which distort the observed CMB temperature 
and polarization fields and so contaminate the B-modes. We identify 7 types of effects, described by 
11 distortion fields, and show their association with known instrumental systematics such as common 
mode and differential gain fluctuations, line cross-coupling, pointing errors, and differential polarized 
beam effects. Because of aliasing from the small-scale structure in the CMB, even uncorrelated 
fluctuations in these effects can affect the large-scale B modes relevant to gravitational waves. 
Many of these problems are greatly reduced by having an instrumental beam that resolves the 
primary anisotropies (FWHM <C 10'). To reach the ultimate goal of an inflationary energy scale 
of 3 x 10 15 GeV, polarization distortion fluctuations must be controlled at the 10 -2 — 10 -3 level 
and temperature leakage to the 10 -4 — 10 -3 level depending on effect. For example pointing errors 
must be controlled to 1.5" rms for arcminute scale beams or a percent of the Gaussian beam width 
for larger beams; low spatial frequency differential gain fluctuations or line cross-coupling must be 
eliminated at the level of 10 -4 rms. 



The recently detected polarization of the cosmic mi- 
crowave background M holds subtle imprints in its pat- 
tern that potentially can reveal the physics of the infla- 
tionary epoch ^ || and provide a new handle on the 
dark matter and energy in the universe [§], §| §]■ This 
curl pattern, the so-called -B-modes, lies at least an or- 
der of magnitude down in amplitude compared with the 
detected main polarization level, which itself is an or- 
der of magnitude lower than the temperature anisotropy. 
Clearly their detection represents a substantial experi- 
mental challenge. 

Beyond raw sensitivity requirements for instruments 
0, much attention has already been given in the lit- 
erature to two aspects of this challenge: astrophysical 
foregrounds (e.g. [S J| |l0)) and the survey mask and 
pixelization (e.g. |l2|). The general requirements im- 
posed on experiments are clear: multiple frequency chan- 
nels and large, contiguous, finely pixelized areas of sky. 
The requirements on other instrumental properties has 
received less attention, in part due to the lack of a com- 
mon language to express their effect on -B-modes. Such 
a language must be expressed in the map, not purely in- 
strument, domain since B- modes reflect a spatial pattern 
of polarization, not its state. In this paper, we seek to 
provide such a connective language and conduct an ex- 
ploratory study on the impact of these systematic effects 
on the science of B-modes. 

We divide polarization effects into two categories: 
those which are associated with transfer between polar- 
ization states of the incoming radiation, mainly induced 
by the detector system (§|), and those which are associ- 
ated with the anisotropy of CMB polarization and tem- 
perature, mainly induced by the finite resolution or beam 
of the telescope (§[H]). We evaluate their effect on S-mode 
science in § [II and on polarization statistics in general in 



the Appendix. 



I. POLARIZATION TRANSFER 

We begin by reviewing the standard transfer ma trix 
formalism for describing polarization detectors in § I A 
and illustrate it s us e in describing the errors in simple 
polarimeters in §IB. The translation to the map domain 
is discussed in §1 C 



A. Description 

The polarization state of the radiation is described by 
the intensity matrix {EiE^j where E is the electric field 
vector and the brackets denote time averaging. As a 
hcrmitian matrix, it can be decomposed into the Pauli 
basis 



P = C(EE t ) 

= ei + Q0-3 



Ucn + Vcr 2 , 



(1) 



where we have chosen the constant of proportionality so 
that the Stokes parameters (<d,Q,U,V) have units of tem- 
perature under the assumption of a blackbody spectrum. 
Note that the Stokes parameters are recovered from the 
matrix as (tr[IP]/2,. . . ,tr[<T2P]/2). We will assume that 
V — on the sky. 

The instrumental response to the radiation modifies 
the incoming state before detection and is generally de- 
scribed by a transfer or Jones matrix J (e.g. jt3), see 0] 
for an introduction in the CMB context), where 



E n 



JE, 



(2) 



The polarization matrix is then transformed as 

Pout = JPin*J^ • 



(3) 



With an estimate of the transfer matrix of the instrumen- 
tal response J, the incoming radiation can be recovered 
as 



J^PoutCjt)" 1 

(J- 1 J)p in (J- 1 J)t 



(4) 



The errors in the transfer matrix determination will then 
mix the determined Stokes parameters according to the 
general transformation rule (^|) with a new transfer ma- 
trix 
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where we parameterized the components with a set of 4, 
possibly complex, numbers (a c , 71 c , 72c, w c ). 

Now let us evaluate the error in the Stokes parameters 
to first order in the real and imaginary parts of the error 
parameters (e.g. Re(a c ) = a, Im(a c ) = 



5{Q ± iU) = (Q± iU) - (Q± iU) 

= (a ± i2u){Q ± iU) + (71 ± i l2 )& ■ 



(6) 



The main effects are a miscalibration of the polarization 
amplitude described by a, a rotation of the orientation by 
an angle u>, and a "shearing" of the temperature signal 
into polarization described by (71, 72), which we will call 
monopole leakage for reasons that will be clear below. 
Note that the imaginary pieces cancel to first order and 
do not appear in Eqn. (|6[). Furthermore terms that cou- 
ple the pair (Q + iU, Q — iU) arise only at second order, 
but we shall see that this need only be true if the Stokes 
parameters are measured through the same transfer sys- 
tem. Note that in the CMB context, monopole leakage 
from (71, 72) is particularly dangerous since the isotropic 
signal is a factor of 10 6 or more larger than the expected 
polarization. 



B. Instrumental Correspondence 

Let us consider a few simple polarimeters that directly 
measure the polarized signal from a single spot on the sky 
(see Fig. [I]; for state of the art techniques, see e.g. [pi). 
Here the incoming radiation is split into two, ideally or- 
thogonal, components, E± and E 2 (using, for example 
an ortho-mode transducer or OMT). These components 
are possibly amplified and coupled into a detector that 
measures the polarization either by differencing or by cor- 
relation. Ideally the transfer matrix of the components 
that split and couple the radiation into the detector is 
proportional to the identity matrix: J oc I. In reality 
it contains systematic errors so that J J cx J. Let us 
parameterize these errors as |fq] 
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FIG. 1: Block diagram for simple polarimeters. The ortho- 
mode transducer (OMT) separates two orthogonal linear po- 
larization states with a leakage between the two characterized 
by (£1,62). After amplification with gain fluctuations (51,32) 
the polarization state is detected by one or more of the fol- 
lowing techniques: differencing the lines to produce Q, cor- 
relating the lines to produce U, correlating the lines with a 
phase shift — -k /2 to produce V . The roles of Q and V may 
be interchanged by placing a quarter-wave plate at the front 
end. 



where g-yp are fluctuations in the gains (or, more gen- 
erally, coupling efficiencies) of the two lines, a is the 
phase difference between the lines, ei 2 express the non- 
orthogonality or cross-coupling between the lines, and 
0i,2 are the phases of these couplings. 

First consider the simple differencing of the time av- 
eraged intensity in the two lines (EyE*) — (EzE^)- This 
forms an estimate of Q in a coordinate system attached 
to the instrument (e.g. Jf^]). Bolometer systems can be 
modeled with this set up, although the means of separat- 
ing the two polarization states and the exact meaning of 
the parameters (gi, 52) may differ (see [|l7| for polariza- 
tion sensitive bolometers). Under the assumption that 
51,2, £1,2, a <C 1, 



SQ = (5l + 92)Q - (f2 COS 02 - £l COS(j>i)U 
+ {gi ~ 92)0 , 



(8) 



so common-mode gain fluctuations act as a normaliza- 
tion error a — (gi + 92) on Q, the cross-couplings act 
as a rotation uj = (e2COS02 — eiCOS0i)/2 and differen- 
tial gain fluctuations leak temperature into polarization 

7i = (5i -£2)- 

Now consider a simple correlation polarimeter where 
the signal in the two lines are correlated as (EyE^) which 
forms an estimate of U in the instrument basis (e.g. p5|). 
Then the errors in the determination become 

SU = (gi + g 2 )U + (e 2 cos0 2 - ti cos0i)Q 

+ (eicos^i +£2cos0 2 )9, (9) 

so again a = (gi + g 2 ), to = (e 2 cos 02 — ei cos0i)/2 but 
leakage from into U is given by 72 = (ei cos 0i + 
e2cos02)- Instead of differential gain fluctuations, the 
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cross-coupling between the lines is responsible for the 
monopole leakage in a correlation system. 

Notice that under the assumption of a <C 1 and van- 
ishing intrinsic V, the phase error a does not appear to 
first order. It is instructive to consider the case where 
a is large, say a — n/2. In this case the correlation 
polarimeter actually measures V not U (see Fig. [l]). In 
general, the phase error a rotates U into V. 

A complex correlation polarimeter actually takes ad- 
vantage of the (U,V) rotation to measure (Q,U) simulta- 
neously (e.g. [18|]). Here, circular polarization states are 
coupled into the lines using, for example a quarter wave 
plate before the OMT. This effectively converts Q into V 
in the instrument basis. The Jones matrix of the quarter 
wave plate is 

\/2 V sin 20 cos 29 — i J ' v ' 

where 9 gives the orientation of the plate with respect to 
the OMT (ideally = 7r/4). After amplification, the sig- 
nal can be coupled into two different correlators, which 
include different phase shifts between the lines. These ad- 
ditional phase shifts can be represented with the transfer 
function: 



' phase 



(0) = 



1 




(11) 



For one correlator <j> is set to zero, yielding an estimate 
of U, while the other correlator has <f> = 7r/2, providing 
an estimate of Q. 

Now let us consider the effect of certain imperfections. 
Consider the actual transfer matrices of the two correla- 
tions to be 

J[/ = Jii„e(5i,S , 2,ei,e2)Ji/4(7r/4 + /3) , (12) 

Jq = Jphasc(7r/2 + ^)Jlino(5l> 52, ei, £2)Ji/4(7r/4 + (3) , 

and the assumed transfer matrices to be 

Ju = J 1/4 (7r/4) , (13) 
Jq = Jphasc(7r/2)Ji/ 4 (7r/4) , 

where the line matrix is taken from Eqn. (Q) with the 
phase factors set to zero for simplicity. Then the errors 
become 

6{Q ± iU) = [( 9l +g 2 )± 2if3] (Q ± iU) 

+ ipU+(e 1 +e 2 )Q. (14) 

The new feature in this system that was not present in the 
simple correlation polarimeter is an asymmetry between 
Q and U which is first order in the phase error ip. More 
generally, a technique that simultaneously measures Q 
and U may have separate transfer properties (calibration, 
rotation, etc.) that appear as a coupling of opposite spin 
states Q + iU and Q — iU. We will call such effects spin 
flip terms. 



The Jones matrix formalism can be applied to more 
complicated polarimeters, such as inteferometers, or 
other systematics such as the finite emissivity of the dish 
polarimeters Jl9| . In general, the systematic errors in 
the detector system will lead to calibration errors, rota- 
tion of linear polarization, leakage of temperature into 
polarization, and coupling between the two spin states 
Q±iU. 



C. Map distortions 

Errors in the polarization sensitivity of the detector 
system that vary with time will translate into errors in 
the polarization sky maps that vary with position. Map 
making generally proceeds by modeling the time ordered 
data as a vector of numbers d (e.g. EQ)) 



d = As + n , 



(15) 



where n is the instrumental noise and s is the model of the 
signal, say {Q(h 1 ),U(h 1 )...,Q(n rip ),U{h np )} for a sky 
map with n p pixels. Here A is the pointing matrix and 
in its simplest incarnation just encodes the sky pixel at 
which the instrument is pointed at the given time. More 
generally the pointing matrix also encodes the beam and 
the chopping strategy where different pointings arc differ- 
enced to remove systematic offsets. The additional com- 
plication for polarization is that the pointing matrix also 
has to encode the orientation of the instrument to trans- 
form (Q,U) in the instrument basis to the fixed sky. This 
is an advantage since systematic errors like the monopole 
leakage (71, 72) are fixed to the instrument basis and not 
the sky. 

Given the statistical properties of the noise N = (nn*), 
the minimum variance map reconstruction is 



s = [A'NA^A'Nd. 



(16) 



This weighting of the data vector then also describes the 
transformation of the instrumental systematic errors to 
errors in the map. In the simplest case of white detec- 
tor noise, fixed instrument orientation, simultaneous Q 
and U detection and no chopping, the weighting sim- 
ply averages the n p separate pointings for each pixel. If 
the systematic fields, e.g. the calibration error a(t) were 
uncorrelated in time, they would remain so in the map 
but with a variance that is reduced by n p . Low frequency 
temporal correlations in the systematics will produce cor- 
related noise in the map. This is generally controlled by 
spatially cross- linking the scans |2l|, HJ. A noise power 
of the 1/ / form will typically lead to spatial correlations 
between white and l/l [f23f , wher e I is the angular fre- 
quency or multipolc moment (sec §111 A). Note that even 
a l/l spectrum gets most of its variance from high / and 
so contamination at the pixel or beam scale will be of 



particular interest in §111 



Since the translation between the temporal and map 
domain is conceptually straightforward but highly de- 
pendent on the scanning strategy, we parameterize the 
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systematic errors directly in the map 

S[Q ± iU](n) = [a± i2w](h)[Q ± iU](n) 

+ [f 1 ±if 2 ](n)[QTiUm 
+ [71 ±i 72 ](n)e(n). 



(17) 



These correspond to calibration and rotation, spin-flip 
coupling and monopole leakage errors as they appear in 
the map. 



II. LOCAL CONTAMINATION 

In the previous section, we dealt with polarization 
transfer in a single, perfectly known, direction on the 
sky. An experiment necessarily has finite resolution and 
thus there is an additional class of contamination associ- 
ated with the resolution or beam of the experiment. We 
will consider here contamination from a local coupling 
between the Stokes parameters which models low order 
anisotropy in the polarized beams. 



A. Description 

Let us consider the polarization fields to be mixed lo- 
cally 

5[Q ± iU}(&; a) = crp(n) • V[Q ± iU](n; a) (18) 
+ a[d 1 ±id 2 ]{h)[d 1 ±id 2 ]Q(n; a) 
+ <j 2 q{n)[d 1 ±id 2 } 2 e{h;a), 

where the fields are smoothed over the average beam of 
the experiment, here denoted by a, the Gaussian width. 
Therefore the fields p, d and q represent sensitivity to 
structure in the fields on the scale of the beam. Wc 
will call these the pointing error, dipole leakage and 
quadrupole leakage respectively. 

We truncate the local coupling at the dipole level for 
the polarization and the quadrupole level for the tem- 
perature. These terms have a direct correspondence to 
known systematics as we shall see. More generally, the 
form of these couplings is dictated by the properties of 
the polarization field under rotation and can be general- 
ized to higher order in a straightforward manner. 



B. Instrumental Correspondence 

Local couplings are primarily due to imperfections in 
the beams. Even a perfectly on-axis, azimuthally sym- 
metric telescope will not in general produce a completely 
azimuthally symmetric and perfectly polarized beam. 
In general the beam has a finite ellipticity along the 
axis of polarization, and there is a "cross-polar" beam 
which couples to the "wrong" polarization state (ad- 
ditional asymmetries may appear in off-axis telescopes 



p5[ ). Both these imperfections arise because the sur- 
face normal to the optics has different orientations with 
respect to polarization axis depending on where the in- 
cident radiation strikes the telescope |l3| , . 

As an illustrative model of these effects, consider the 
case where the receiver on the telescope is a simple dif- 
ferencing polarimeter. Here the effects of the cross-polar 
beams are of second order, and will be ignored. Radi- 
ation from the sky is then coupled into one line of the 
detector through a perfectly polarized beam: 



B(n;b,e) = 



1 



2ncr 2 (l-e 2 ) 

, ( n 2 



exp 



1 

2^ 



(ni - h) 2 
(1 + e) 2 



b 2 f 



{1-ef 



(19) 



where b is the offset between the beam center and the 
desired direction on the sky, a is the mean beamwidth, 
and e is the ellipticity These parameters are differ- 
ent for the different polarizations, and the difference in 
the beams enters into the Q measurement 



B(n;b ,e ) - B(n;h b ,e b ) . 



(20) 



To first order in the sums and differences of the elliptici- 
ties and pointing errors 



(jp = (b„ + b 6 )/2 , 
ah d = (b Q - b 6 )/2 , 

e s = (e„ + e b )/2 , 
q = (e a - e b )/2 , 



(21) 



we obtain 



Q{h-a)= J dn'S(n'){Q(n + n' + crp) 



b d n' 



a a 

Q(n; a) + crp • VQ(n; a) + ab d • V9(n; a) 
+ a 2 q[d 2 1 -d 2 2 Mh;a) 



e(n + n')}, 
r) 

(22) 



where the average beam B(h) = B(n;0,0) and we drop 
second derivative terms in Q. A difference in the mean 
beamwidth of the two beams has the same form as a con- 
tribution to the monopole leakage except that the filter 
for the temperature field is the beam difference not beam 
and is not simply a low pass filter. 

A pointing offset in both beams becomes a gradient 
coupling in polarization and a differential beam elliptic- 
ity or "squash" becomes a coupling to the temperature 
quadrupole. These have a clear correspondence to the 
local contamination model of Eqn. (p"8|). A differential 
pointing offset, or "squint" translates into coupling to 
the temperature dipole. The exact correspondence to the 
model is not precise since the leakage does not truly be- 
have as a false polarization. For example under rotation 
of the instrument by 7T, the false Q reverses sign. The 
model of Eqn. (Il8|) with d referenced to the sky (not the 
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instrument) does transform as polarization and so should 
be viewed as the residual dipole sensitivity after correc- 
tion. The quadrupolar coupling is particularly dangerous 
since it behaves precisely as a polarization and cannot be 
removed through rotation of the instrument. For exam- 
ple, even a circularly symmetric temperature hot spot 
becomes a radial pattern of polarization through its lo- 
cal quadrupole moment. 

These leakage terms also appear if the receiver is a 
correlation polarimeter. In this case, the leakage from 
temperature to polarization is due to the amplitude and 
shape of the cross-polar beam instead of asymmetries in 
the main beam However, because both of these im- 
perfections have a common origin in the variations in the 
boundary conditions at optical surfaces, it turns out that 
a differencing system or a correlation polarimeter will ob- 
tain similar leakage terms due to the local effects (after 
accounting for the orientation of the receiver with respect 
to the optics) |24|, |2?|| . Interferometric polarimeters have 
related effects on the scale of the primary beam although 
some of the implications for B-modes will differ since the 
measured modes are below the beam scale [p8| . 

If stable, these effects can be removed given a beam 
measurement and the true temperature field on the 
sky using the formalism of anisotropic polarized beams 
p9j |30f . Moreover, as the circularly symmetric hot spot 
example implies, a stable quadrupole le akage produces 
no .B-mode in the polarization map (see § 111 B ). It is the 
instability in these effects or errors in the subtraction 
that appear as errors in the map. 



III. B-MODE CONTAMINATION 

We study here the implications of polarization trans- 
fer and local contamination on the i?-modes of the po- 
larization. In § III A , we give the harmonic representa- 
tionof the polarization and contamination fields. In 



III B, we compute the contamination to the B power 



spectrum from polarization distortion and temperature 
leakage. We explore the implications of these effects in 



IIIC 



A. Field Representation 

The polarization and contamination fields may in gen- 
eral be decomposed into harmonics appropriate to their 
properties under rotation or spin. For small sections of 
the sky, these harmonics are simply plane waves jn], ||| ; 
in the Appendix we treat the general all-sky case. We 
will follow the convention that a complex field S of spin 
±s is decomposed as 



d 2 i 

(2^ 



[S a ±iS b }(\) 



(23) 



where cos</>; = l x /l. The complex polarization Q ± ill 
is a spin ±2 field and we will follow the conventional 




FIG. 2: Scalar CMB power spectra in temperature (00) and 
_E-mode polarization (EE) compared with _B-mode polariza- 
tion due to gravitational lensing and gravitational waves at 
the maximum allowable 2.6 x 10 16 GeV |36| and minimum de- 
tectable 3.2 x 10 15 GeV level EtIi. The ACDM model shown 



has parameters given in III A 



nomenclature that its harmonics are named E ± iB. 
This property requires the calibration a, rotation lu, 
and quadrupole leakage to be spin-0 fields, the point- 
ing pi ± ip2 and dipole leakage d± ± id,2 to be ±1 fields, 
the monopole leakage 71 ± 172 to be ±2 fields, and the 
spin-flip fi ± if 2 to be ±4 fields. For spin ±1 fields S a is 
the divergence-free part and St is the curl- free part. 

Under the assumption of statistical isotropy of the 
fields, their two point correlations are defined by their 
(cross) power spectra 



<S(1)*5'(1')) = (2n) 2 S(l - l')Cj 



ss' 



(24) 



where S, S' are any of the fields. In particular, the CMB 
polarization is described by Cf E and C BB and CMB 
temperature by Cf . Note that Cf B — for scalar fluc- 
tuations in linear theory. For definiteness, let us take as 
a fiducial model: a baryon density of VLbh? = 0.02, cold 
dark matter density of fl c h 2 — 0.128, a cosmological con- 
stant of Ha = 0.65, reionization optical depth r = 0.05, 
an initial amplitude of comoving curvature fluctuations of 
or erg = 0.92), and a scalar spectral 



4.79 xl0~ 5 ( 33 



index of n = I in a spatially flat universe. Power spectra 
for the fiducial model are shown in Fig. |^. It is the large 
range in expected signals that make the contamination 
problem for BB so problematic. 

We will calculate the contamination to the -B-mode po- 
larization power spectrum assuming no intrinsic i?-modes 
and generally will plot 



AB 



1(1 



2- 



l \ C BB 



1/2 



(25) 



in units of fiK. The general case is given in the Appendix. 
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(a) Calibration 



-lensing 



(b) Monopole Leakage 



-lensing 




FIG. 3: Coherence dependence of B-mode contamination (a) for calibration a with rms A a = 10 



7a, 7t with A la — A lb = 10 added in quadrature. The beam scale is FWHM 



!ln2) 



1/2, 



2 (b) for monopole-leakage 
1' to remove beam effects 



and the FWHM coherence (81n2) 1//2 a is stepped from 256' to 4' in factors of 2. Other effects follow the trend of calibration 
errors not monopole leakage. For a coherence large compared with the CMB acoustic peaks, B contamination picks up their 
underlying structure. Here and in the following figures, the gravitational lensing and minimum detectable gravitational wave 
(Ei = 3.2 x 10 15 GeV) B-modes are shown for reference (thick shaded lines). The scaling with Ei of the peak in the B-mode 
spectrum is shown on the right hand axis. 



Although the distortion fields need not be statistically 
isotropic, for illustrative purposes we will take contami- 
nation fields with power spectra of the form 



Cj= b oc exp(-Z(Z + l)c4), 



(26) 



i.e. white noise above some coherence scale as- The 
normalization constant is set so that 



A 2 
A s 



(27T) 



ss 

2^1 ' 



(27) 



The set {As, as) then characterizes the rms and coher- 
ence of the contamination field. 



B. JS-modes 

The changes to the -B-mode harmonics due to the cal- 
ibration, rotation, spin-flip and pointing take the form 



5B{\) 



d 2 h 



S(h)E(h)W s (h,h) 



(28) 



with I2 = 1 — li and 



W a = sm[2{fa 2 - fa)} > 
W u = 2cos[2{fa 2 -fa)}, 
W Pa = er(l 2 x l x ) ■ zsin[2(^ 2 - fa)] , 

W Pb =a{h-h)sm[2{fa 2 -fa)], 

W fa = sm[2{2<f> h - fa 2 -fa)] , 
W fb = cos[2(20 ;i - fa 2 - fa)} , 



(29) 



for the various effects. Here li = Zilj.. These relations 
imply contamination to the BB power spectrum of 

*Cf B = E/ &<%*Cg*{*)WSWs, , (30) 



where 



C^(a)=C^eM-Kl + l» 



(31) 



is the EE power spectrum smoothed over the average 
beam. 

Similarly the change due to temperature leakage can 
be described by 



SB(l) 



J ^S(h)e(h)W S (h,h), (32) 



with 



W la = sin[ 2 (<f>h -<h)] , 
W lb = cos[2(0 /l - fa)} , 
W da = -{her) cos[(t) h + fa 2 - 2 fa] , 
W db = {l 2 a) sin[0 ;i + fa 2 - 2fa] , 
W q = ~{l 2 <j) 2 sm[2{fa 2 -fa)}, (33) 

leading to 

6C " B = E / ^C?fC?°{a)W* s W s . , (34) 

SS' ^ ' 

for the power spectrum contamination ]34| . 

A few limiting cases are worth noting before proceed- 
ing to specific examples. If l\ ^> I as is the case for 
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(a) Pointing 



g-lensing 



(b) Quadrupole Leakage g-lensing 




FIG. 4: Beam dependence of B-mode contamination for (a) pointing with an rms A PtL = A Pb = 10 2 (in units of the Gaussian 
beam width) added in quadrature (b) quadrupole leakage with an rms A q = 0.002 (in units of differential beam ellipticity) . 
The coherence a is set to max(o", 10'/(8 In 2) 1 ' 2 ) and the beam is stepped from 128' to 2' in factors of 2. 



power in the contamination field at much smaller scales 
than the I of interest, l\ rs li and fa w — fa. The geo- 
metric factors in Eqn. ( p9| ) and (]3^) cause all effects to 
efficiently produce S-modes except the pointing curl p a 
where the cross product vanishes. In the opposite limit 
h <C I then 1% w I and fa ~ 0. Here the calibration a, 
pointing terms, and quadrupole leakage are geometrically 
suppressed. The reason is clear from the nature of the 
effects: a uniform distortion in any of these quantities 
does not produce a B-mode. 



C. Scientific Impact 

Cosmological B-modes come from two main sources: 
gravitational waves, also known as tensor perturbations, 
|^|, |j| and gravitational lensing of polarization by the 
large-scale structure of the universe . Aside from small 
but interesting effects due to the dark energy, reioniza- 
tion and massive neutrinos, the gravitational lensing B- 
modes can be predicted given parameters extracted from 
the CMB temperature spectrum. The gravitational lens- 
ing prediction in the fiducial model is shown in Figs. |||5| 
as the shaded top line. 

Under slow-roll inflation, the initial amplitude of the 
gravitational wave spectrum is parameterized by the en- 
ergy scale of inflation Ei and its spectrum is nearly scale 
invariant. It predicts a i?-mode power spectrum ampli- 
tude with a peak at I ~ 90 of [p5[ 



A£ pcak = 0.024 



10 16 GeV 



(35) 



Under reasonable cosmological assumptions, the CMB 
temperature anisotropies constrain the energy scale to be 
Ei < 2.6 x 10 16 GeV (36). If the energy scale is less than 
Ei < 3.2 x 10 15 GeV, then even with a direct reconstruc- 
tion of the lensing signal 0, a significant detection of 



the inflationary B-modes cannot be achieved |$7) . These 
two extremes are shown in Fig. ||and are used in Figs. 
||, to mark the range across which the systematic errors 
need to be controlled. We will take the prediction for 
the middle of this range (Ei — 10 16 GeV) as the minimal 
level that a next generation polarization mission must 
reduce errors. For reference, with no systcmatics or fore- 
grounds the Planck satellite ]3q | can in principle achieve 
a la bound of E t = 1.1 x lO^GeV §§. 

Given that the inflationary £?-modes peak at I sa 90, 
one might naively assume that only contamination fields 
with coherence corresponding to degree scales would be 
problematic. However because equations (|2^ ) and (|3^ ) 
represent mode coupling, this expectation is incorrect. 
The problem is that the intrinsic power in the CMB po- 
larization fields as well as the temperature gradient and 
second derivative fields peak on the scale associated with 
the diffusion scale at recombination, now observationally 
determined to be / w 10 3 or 10' by the CBI experiment 
p9[ . In Fig. ||a, we show the effect of a calibration error 
with the same rms A a = 10~ 2 but different coherence 
scales a a . For coherence scales above (81n2) 1 / 2 a a = 10', 
the contamination actually increases as the coherence 
scale decreases. For most effects, the coherence scale that 
gives the maximum total contamination is the larger of 
the beam scale and ~ 10'. The mathematical reason is 
that the mode coupling sets 1 = li + I2 which forms a 
triangle with sides (l,h,h)- For CMB power at 1 2 ^> 1 
contamination power at l\ w l 2 causes most of the leakage 
by forming a flattened triangle. 

The exception is the monopole leakage which takes 
power out of the CMB temperature power spectrum it- 
self, not derivative power spectra which are weighted by 
factors of I. Here the most damaging coherence scale is 
associated with the first peak in the CMB at I w 200 
(see Fig. ||) which is dangerously close to the I ~ 100 
scale of interest for gravitational waves. Fig. ||b illus- 
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Cs (f) ps W 


C S (2°) p s (2°) 


Calibration a 


0.060 


-0.3 


0.049 


0.0 


Rotation w 


0.015 


-0.3 


0.011 


0.0 


Pointinp" T> r . 


0.75 


-1.3 


0.53 


-1.0 


Pointing p/j 


0.098 


-0.7 


0.57 


-1.0 


Flin f„ 

1 11 P Jo, 


0.061 


-0.3 


0.046 


0.0 


Flip f b 


0.059 


-0.3 


0.045 


0.0 


Monopole f a 


0.0023 


-0.9 


0.0006 


0.0 


Monopole 75 


0.0019 


-0.9 


0.0005 


0.0 


Dipole d a 


0.0077 


-1.3 


0.0053 


-1.0 


Dipole d;, 


0.0077 


-1.3 


0.0056 


-1.0 


Quadrupolc q 


0.0124 


-1.5 


0.0394 


-2.0 



TABLE I: Scaling parameters for contamination effects with 
a coherence of the beam scale a and 2°. Cs represents the 
minimum rms required to not exceed a signal at Ei = 10 16 
GeV; the ultimate limit of 3.2 x 10 GeV would require an 
order of magnitude smaller rms. 



trates this problem and shows that degree scale fluctua- 
tions in monopole leakage from low frequency noise must 
be controlled to substantially better than 1CP 3 rms for 
Ei < 10 16 GeV. 

Pointing, dipole and quadrupole leakage errors are ex- 
pressed in terms of fractions of the beam and hence can 
depend strongly on the beam scale. The contamination 
at I ss 90 from pointing errors of a fixed fraction of 
the beam holds roughly constant for beam sizes above a 
FWHM « 10'. At this point most of the structure in the 
underlying CMB fields become resolved and the contam- 
ination depends on the absolute pointing error relative 
to the CMB 10' coherence (see Fig. |J). Pointing prob- 
lems must be constrained to better than the larger of 
10 _1 of the Gaussian beamwidth or 15" absolute rms for 
Ei < 10 16 GeV and must reach 10 -2 of the beamwidth or 
1.5" absolute to be safely irrelevant. 

The quadrupole leakage provides a more extreme ex- 
ample. Contamination for a fixed rms differential ellip- 
ticity strongly increases with increasing beam and so a 
beam < 10' FWHM greatly reduces the contamination. 
The dipole leakage lies in between these two cases in sen- 
sitivity to the beam scale. In Fig. || we show all of the 
effects, for a choice of beam and coherence of FWHM 
= 10' and an rms of 10~ 2 for polarization distortions 
and 10~ 3 for temperature leakage. 

It is useful to have an approximate scaling for the rms 
amplitude of the systematic needed to make the contam- 
ination on the same level (at I = 90) as a given target 
inflationary energy scale. Let us approximate the rms as 
a power law in the FWHM of the beam 

As=Cs {W*Gev) ( ^0^ } ■ (36) 

In Table 1, we give the coefficients Cs and ps for 
two choices of the coherence scale: as — o~ and 
(81n2) 1 / 2 a5 = 2°. The beam dependence is calculated 
locally around 10' and should not be used to extrapolate 
results far from this. 



IV. DISCUSSION 

We have provided a fairly general description of the 
phenomenology of systematic errors that can occur in po- 
larization maps, their correspondence with known classes 
of instrumental problems, and their impact on the sci- 
ence of .B-modes. Instability in the systematic effects 
or errors in their removal lead to residual contamination 
in the polarization maps that are parameterized by 7 
fields, 4 of which have two components each, for a total 
of 11 distortion parameters per position on the sky or 
multipolc moment. These errors are associated with cal- 
ibration, rotation, pointing (2), spin flip (2), monopole 
leakage (2), dipole leakage (2) and quadrupole leakage. 
The three temperature leakage effects are named for the 
type of temperature fluctuation across the beam scale 
that they respond to and are especially dangerous due to 
the extremely low level of polarization expected in the B- 
modes. Monopole leakage generally arises in the receiver; 
dipole and quadrupole leakage are associated with asym- 
metries in the beam. 

We have illustrated these problems by modelling the 
fluctuations in these contamination fields with an rms 
amplitude and coherence. In general, it is not sufficient 
to control the fluctuations in the field on the degree scales 
of interest for gravitational wave B-modes. Because all of 
these effects transfer power from the CMB fluctuations 
themselves, the most dangerous fluctuations are those 
that are on the same scale as most of the power in the 
CMB fields. For all but the monopole leakage effect, 
which can draw power out of the first acoustic peak, the 
underlying power lies at the diffusion damping scale of 
I ~ 10 3 or ~ 10'. Unless the beam resolves this scale, 
even uncorrelated white noise fluctuations in the fields 
can substantially contaminate low multipoles in B. 

The interplay between the beam scale and 10' coher- 
ence scale of the CMB fields plays an especially important 
role in pointing, dipole leakage and quadrupole leakage. 
These problems couple local derivatives of the CMB fields 
into false polarization signals. They can largely be elim- 
inated if the beam is sufficiently small so that the CMB 
fields are smooth across the beam scale. Small beams are 
also desirable for constructing weak lensing mass maps 
from the -B-mode polarization Q . 

Based on the systematic errors of the current genera- 
tion of experiments, these problems should be challeng- 
ing but not insurmountable. The DASI instrument had 
percent level monopole leakage which was stable at the 
fractional percent level and quadrupole leakage also at 
the percent level which was highly stable. It also had 
rotational uncertainties at the percent level ^8|. The 
PIQUE instrument had a monople leakage at under the 
percent level and a dipole leakage at less than 2.5% |fto| . 
The polarization sensitive bolometers for the upcoming 
Boomerang experiment |l7| and planned for Planck have 
monopole leakage at the percent level but are claimed to 
be very stable. 

This exploratory study of polarization effects should 
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(a) Polarization Distortion g-lensing 




(b) Temperature Leakage g-lensing 

monopole (y a , y 

dipole (d a , djj) 

quadrupole (q) 



FIG. 5: All effects for a beam and coherence of FWHM = (81n2) 1/2 cr = 10'. (a) Polarization distortion for an rms of A = 10~ 2 
from calibration a, rotation u (0.6° rms), pointing (p a ,Pt) (2.5" rms) , and spin flip (fa, ft)- (b) Temperature leakage for an 
rms of A = 10~~ 3 from monopole (-y a , 75), dipole (d a ,d ) and quadrupole (q) terms. The "b" component of each effect is shown 
with dashed lines. 



help to provide some rough guidance on the long road 
ahead toward the ultimate goal of detecting the gravita- 
tional waves from inflation. 
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APPENDIX A: GENERAL TREATMENT 



where X + = E, X = iB and S G a,aj,p a ,pt, f a , /& and 



The flat sky expansion in Eqn. fl23|) may be generalized 
by decomposing the fields as [ l| 

[Si ± iS 2 ](n) = (i) s J2i S a ± iS b ]lm ±sY lm (h) , 

In i 



where s Y/ m is the spin-s spherical harmonic |41j| . 

The corrections to the E and B harmonics of the po- 
larization from the distortion field S can be generally 
expressed as 



sxt = (-ir E E 



(2J + l)(2Ji + l)(2l 2 + l) 



-in 



(Al) 



1 -m mi m 2 ' Sl ^ L s 



hm 2 S / 2 iri2 ' ' 



X hm, ) 



a,p b 

4. 



i[i±(-i)'^+«»], 



(A2) 



selects out even and odd sums of the Vs. The factor 
<Zg = 1 for S G a,ui,p a ,Pb and = ±1 for S G f a , fb and 
adjusts the relative sign. 

The specific linear source terms are 



_ _i _ I I hh 



L n „ — —— 



V(l2 + 2)(l 2 -l) 



- V(*2-2)(Z 2 +3) 

■ T ( I h h 

& = -2 4-2 



I h h 
-2 1 1 

I h h 
-2 -1 3 

I h h 
-2 1 1 

I h h 
-2 -1 3 



(A3) 



Under the assumption of statistical isotropy of the dis- 
tortion fields, the perturbation to the power spectra are 
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given by 



for S e "ta,7b,d a ,d b ,q 



(2; 1 + i)(2; 2 + i) c g S , 

' \/47r ! ^ 4tt (l 

v hl 2 SS> 

x 44^44 + C* B e~ s e- s ,]L* s L s> , 

BB a 00 ,^BB , ( 2 ^1 + 1 )( 2 ^2 + l^SS' 



47T 



x o c a 



2w, 



'00 / r iEE nBB\ 



§C EB = ^L( C EE _ C BB- ) 
V47T 



(A4) 



E 



(2/ 1 + l)(2i 2 + l) ^ sg < 
4n 6(1 



;i/ 2 ss' 

x 4«s'[Qf B 4e S ' + C? 2 B e- s e+){-i)L* s L s , , 

where we have allowed for the possibility of a monopolc 
term in the calibration a and rotation u>. Such terms 
often arise from second order terms in an expansion and 
comes about through the variance of a distortion field 
across the sky. They must be kept since the net change 
to the power spectrum is itself second order and often 
cancel the linear effects. For example, second order terms 
in the pointing errors appear as a monopolc calibration 
error. Since these terms do not transfer power in EE to 
BB, they are not relevant for the discussion in the main 
paper. 

Temperature leakage terms may similarly be described 
in their effect on E and iB 



"lb, da 2 



1 [i± (_!)«-*+«,] ) 



[1T(-1) 



i+h+h 



(A6) 



and 



t» ~ to ~ I -2 2 



L db = iL da = ^l 2 (h + l)a ( J 2 \ \ ) , (A7) 



'(i + 2)! , f I hi 



q y (1-2)!" V" 2 2 J ' 
The perturbation to the power spectra are given by 

5C? E = J2 {2ll + 1){2l2 + 1} Cf s 'cf e 



ihSS' 

e%e%,L* s L s > , 



An 



scr= e {2h+1 f l2+l) c ™' c ® @ (as) 

hl 2 SS' 

e s e s ,L *s L s> , 
§(jEB _ (2/i + 1)(2Z 2 + 1) ss > ee 

hhSS' 

e s e s'(- i ) L *s L s' ■ 



s*L = E E 



(2Z + l)(2/i + 1)(2Z 2 + 1) 



47T 



(A5) 



y 1 -m mj m 2 J A, '"» e, ™ e S is ' 



This completes the general description of the polarization 
contamination from the class of map distortions consid- 
ered. 
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